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©IMETHODE POUR REPRESENTER DANS UN DOMAINE SPATIO-TEMPOREL LA TRAJECTOIRE D'AU MOINS 
^ UN purrs DE FORAGE. 

, 67) 1) Methods pour representer dans un domaine spatio- 
remporel la trajectoire d'au moins un puits de forage. 
2VLa methods est caracterisee en ce qu'elle consiste a 

a) ^laborer un module vitesse-profondeur (4) cal6 au 
puits (1), geologiquement representatif de ladite zone et 
comportant une representation de ladite trajectoire (5) dudit 
puits; 

b) selectionner au moins un point (A) sur ou pr^s de la- 
dite trajectoire, (edit point (A) etant defini par ses coordon- 
nees spatiales et par au moins un couple de donn^es 
pendage-azimut; et 

c) au moins mod6liser dans un domaine somme synth^ti- 
que (9) I'image (a) du point selectionne (A), ladite image (a) 
6tant d'une part d§fmie par des coordonn^es spatio- 
temporelles et par un vecteur pendage-temps et, d'autre 
part, obtenue par un lancer de rayon a deport nul dans le 
modele vitesse-profondeur (4). 
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METHODE POUR REPRESENTER DANS UN DOMAINE SPATIO- 
TEMPOREL LA TRAJECTOIRE D'AU MOINS UN PUTTS DE 

FORAGE 

La presente invention conceme une m6thode pour representer dans 
un domaine spatio-temporel la trajectoire d'au raoins un puits de forage et 
plus particulidrement pour caler dans un bloc de donn^es sismiques 
migr6es-temps la trajectoire d'un puits devi6. 

Le calage d*un puits dans la sismique 2D ou 3D est necessaire 
lorsqu'on souhaite projeter sur une reprdsentation sismique a deux ou trois 
dimensions du milieu dans lequel le puits a ete forfi toutes les informations 
relatives a des donnees diverses qui sont enregistr^es localement dans ou a 
proximity du trou de forage, generalement appel6es logs. II en est ainsi 
des informations relatives aux logs sonique, de densite, de r6sistivite, etc., 
a la profondeur des couches geologiques dudit milieu, aux limites 
biostratigraphiques, aux pendages des r6flecteurs definis notamment par 
les interfaces entre deux couches consecutives, aux azimuts ou meme aux 
images foumies par les diff6rentes techniques de sismique de puits. 

Le calage a done lieu au moins 

a) lorsque dans le traitement des donnees sismiques, le ou les 
puits sont utilises comme r6f6rence dans une d6convolution i phase-z6ro, 
une calibration AVO (Amplitude Versus Offset) ou lorsque les vitesses 
dans le puits sont utUisees en combinaison avec les vitesses de sommation 
pour construire des modeles de migration, 

b) lorsque Tinterpretateur commence une interpretation pour 
connaitre quel horizon sismique correspond a un marqueur ou reflecteur 
geologique donn6, 

c) lorsque les representations profondeur produites dans une 
conversion temps-profondeur, doivent etre calees au puits en termes de 
profondeur et de pendage. 

Le calage pent etre realis6 soit dans le domaine profondeur, soil 
dans le domaine de vitesses d'intervalle, soit encore dans le domaine 
temps. 

Le calage dans le domaine profondeur signifie qu'on souhaite 
disposer d'un modele profondeur ou d'une s6rie de cartes (ou 
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repr6sentations) profondeur qui sont coh6rentes avec la profondeur k la 
position du puits. 

Le calage dans le domaine de vitesse d'intervalle est souhait6 
lorsqu'on ddsire avoir un modele de vitesse qui comprend des vitesses 
d'intervalle coh6rentes avec les vitesses d^intervalle estimees sur les puits 
par un log sonique calibre ou toute autre m6thode. Ceci est generalement 
le cas pour les modeles utilises pour la migration. 

Lorsque le calage est realis6 dans le domaine temps, la trajectoire 
du puits, qui est d6finie en profondeur, doit etre reprdsentee dans ce 
domaine temps. Or, c'est reellement ce calage qui pose un veritable 
probleme. 

En effet, lorsqu'on d6sire comparer une s6rie d'attributs ou de 
valeurs attachees a chaque point en profondeur de la trajectoire du puits et 
relatives par exemple au log sonique, densite, electriques, aux limites de 
couches, aux pendages, aux azimuts, etc., il serait ideal de repr6senter 
i' image sismique en profondeur, c'est-a-dire realiser une migration 
profondeur en 3D, et de faire le calage dans ce domaine "naturel". 
Cependant, cela exige une connaissance precise du champ de vitesse du 
milieu. En pratique, on realise une projection de la trajectoire du puits 
dans le domaine temps-migre (X, Y, T). Pour ce faire, on suppose que la 
migration a deplac^ T information k sa position correcte X, Y et que le 
temps migre est 6gal a un temps de trajet vertical. II s'ensuit que le log 
dans le domaine temps est une version 6tir6e verticalement du log en 
profondeur, la relation profondeur-temps 6tant d6riv6e en general de 
r integration du log sonique calibre par sismosondage (check-shot survey 
en Anglais) ou par les temps des premidres arriv6es d'un PSV (Profil 
Sismique Vertical) par exemple. Malheureusement, retirement vertical et 
le repositionnement par la migration temps en X, Y sont des 
approximations trop larges ou grossieres quand un resultat pr6cis est exig6 
(par exemple pour la delineation d'un reservoir) ou dans les cas d'un 
milieu complexe ou tectonis6. 

Comme on peut le voir, une premiere operation essentielle est 
I'dlaboration d'un modele de vitesses de propagation en profondeur. 

Plusieurs techniques sont utilises ou ont ete proposees. Pamii 
celles qui sont le plus couramment utilis6es, on peut citer celle relative k 
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la conversion veiticale temps-profondeur, celle utilisant la migration le 
long du rayon-image (correction dite de HUBRAL), celle utilisant la 
migration le long du rayon nonnal, celle utilisant la tomographie ou 
incluant T inversion des vitesses de sommation et celle utilisant le proc£d6 
ditde ''coh6rence". 

La conversion verticale temps-profondeur et la migration le long du 
rayon image sont effectuees a partir de reflecteurs migras-temps et 
d' informations sur les vitesses de propagation des ondes sismiques 
extraites soit du puits lui-mame (sismosondages, PSV, etc.) soit de 
donnees sismiques de surface (vitesses de sonmiation ayant subi 
r inversion par application de la formule de DIX). 

La migration le long du rayon normal utilise des r6flecteurs dans le 
domaine de sommation, lesdits reflecteurs 6tant obtenus soit par un point6 
dans le domaine sommation, soit par une d6migration des reflecteurs 
migres-temps. Cette migration utilise le meme type d' informations sur les 
vitesses de propagation des ondes que celui rappeI6 ci-dessus. 

Toutefois ces procddes pour rdaliser une carte ou representation 
profondeur presentent de nombreux inconvenients. En particulier, les 
vitesses de sonmiation sont dans le domaine sommation-temps alors que 
les horizons sont dans le domaine temps-migr6. Or, les deux series de 
donn6es ne sont pas coh6rentes. De plus, les vitesses de sommation 
comprennent des informations sur la vitesse d'intervalle autour du rayon 
normal plutot qu* autour de la verticale. Un autre inconvenimt reside dans 
le fait que les vitesses de sommation sont affectdes par le pendage des 
horizons, les variations lat6rales des vitesses et par d'autres facteurs. 
Certes, on peut remedier k certains de ces inconvenients en appliquant un 
traitement DMO (DIP MOVE OUT en Anglais), qui est mis en oeuvre sur 
les collections en points milieux communs (CMP) avant Tanalyse de 
Vitesse de sommation. 

Ainsi, les proced6s tels que : inversion des temps de trajet, 
inversion de vitesses de sommation, inversion par coherence ou d'autres 
techniques plus 61abor6es faisant appel k la migration profondeur avant 
sommation, peuvent Stre utilis6es pour I'daboration des moddles vitesse- 
profondeur. 
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Pes proc6d6s ont 6t6 con^us pour r6duire les approximations 
constat6es. Cependant, dans leur version classique, le modde derive de la 
Vitesse de sonunation est un module isotrope, et il ne comprend toujours 
pas de moyens pour effectuer un calage du puits, en temps migre. Le 
5 produit final obtenu avec ces proc6d6s est une serie de cartes profondeur 
et des intervalles de vitesse ou fonctions de vitesse [Vo(x,y) + KZ] qui 
sont obtenues seulement a partir de la sismique. 

Lx)rsque ces cartes sont obtenues a partir seulement de donnees 
sismiques de surface, on doit s'assurer qu'il y a coherence entre la 

10 profondeur mesur^e sur les puits et les cartes profondeur. Cela est 
classiquement fait en calibrant des cartes de vitesse d'intervalle par 
Tutilisation des vitesses d'intervalle resultant par exenn{)le des relations 
temps-profondeur mesurees a Templacement du puits. Un autre proced6 
est de caler directement les cartes en profondeur. 

15 Un autre inconv6nient de ces proced6s reside dans le fait que 

r anisotropic du sous-sol est prise en compte de manidre trds approchde 
dans la determination des vitesses. 

La pr6sente invention a pour but de proposer une methode pour 
caler la trajectoire d'un puits de forage d6vie dans le domaine temps- 

20 migre en amaiiorant, en temps reel, le modele de vitesse-profondeur, ce 
qui est Tobjectif recherche. 

La prdsente invention a pour objet une methode pour repr^senter 
dans un domaine spatio-temporel la trajectoire d'au moins un puits devie 
et plus particulierement une methode de calage de ladite trajectoire qui 

25 inclut un maximum d'heterogeneites et Tanisotropie du milieu dans lequel 
le puits devie est fore, laquelle methode est caracterisee en ce qu'elle 
consiste a : 

a) eiaborer un modele vitesse-profondeur caie au puits, 
geologiquement representatif de la zone et comportant une representation 

30 de la trajectoire dudit puits , 

b) seiectionner au moins un point sur ou prds de ladite 
trajectoire, ledit point etant defini par ses coordonnees spatiales et par au 
moins un couple de donnees pendage-azimut ; et 

c) au moins modeiiser dans un domaine somme synthetique 
35 rimage du point seiectionne, ladite image etant d'une part definie par des 
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coordonntes spatio-temporelles et par un vecteur pendage-temps et, 

d'autre part, obtenue par un lancer de rayon k d6port nul dans le moddle 

vitesse-profondeur, 

Selon une autre caract6ristique, la methode comprend en outre 
5 Tetapc d'effectuer une migration en temps de ladite image du domaine 

synth6tique , de maniere a obtenir une repr6sentation de ladite image dans 

un domaine migre-temps. 

Selon une autre caract6ristique, le domaine somme synthdtique 

contenant T image est superpos6 a un bloc somme de donnees sismiques. 
10 Selon une autre caracteristique, le domaine migre- temps contenant 

la representation est superpose a un bloc de donnees sismiques migrees- 

temps. 

Selon une autre caracteristique, le bloc de donn6es sismiques 
migrees-temps est obtenu a partir du bloc somme synthetique en utilisant 
15 un modele de vitesse migration donn6 et en ce que la migration de ladite 
image (a) est effectu6e au moyen dudit modele de vitesse. 

Selon une autre caracteristique, le module de vitesse-profondeur est 
k couches delimitdes par des r6flecteurs. 

Selon une autre caracteristique, le moddle de vitesse-profondeur est 
20 un moddle maill6. 

Selon une autre caracteristique, on selectionne plusieurs points 
dans le moddle vitesse-profondeur et plus particulidrement ies points situ^s 
proches de intersection de la trajectoire avec les r6flecteurs. 

Selon une autre caracteristique, le bloc de donnees sismiques 
25 migrees-temps est obtenu par migration de KIRCHHOFF du bloc de 
donnees sismiques somme et en ce que le domaine migre-temps resulte 
d'une sommation analytique de ladite migration de KIRCHHOFF. 

Selon une autre caracteristique, le bloc de donnees sismiques 
migrees-temps est obtenu au moyen d'une migration par differences fmies 
30 et en ce que chaque image dans le domaine somme synthetique est 
distribuee sur une surface dite surface de diffraction. 

Selon une autre caracteristique, on realise au cours de la mise en 
oeuvre de ladite methode une comparaison entre le domaine somme 
synthetique et le bloc somme enregistre et on determine I'ecart a la 
35 tangence ou ^ la coincidence susceptible d'apparaitre entre chaque surface 
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de diffraction et le r6flecteur correspondant et, lorsqu'il y a un ficart, on 
procMe Si des iterations successives sur le moddle vitesse-profondeur 
jusqu'i cc que ledit 6cart soit acceptable pour consid6rer que ledit modde 
est correctement cal6. 

Selon une autre caract6ristique, la representation en temps-niigrg de 
la surface de diffraction est une surface migree-temps. 

Selon une autre caracteristique, on d6tennine, au cours de la mise 
en oeuvre de ladite methode s'il y a un ecart a la tangence ou a la 
coincidence entre ladite surface migree-temps et T horizon migre-temps 
correspondant et, lorsqu'il y a un ecart, on procede a des iterations 
successives sur le modele vitesse-profondeur jusqu'a ce que ledit 6cart soit 
acceptable pour considerer que ledit modele est correctement cal6. 

Selon une autre caract6ristique, le point de tangence ou de 
coincidence (a ou a*), lorsqu'il existe, permet de determiner le pendage 
reel du reflecteur dans le modele vitesse-profondeur. 

D'autres avantages et caracteristiques apparaitront k la lecture de la 
description d'un mode de realisation de la pr6sente invention, ainsi que 
des dessins annexes sur lesquels la figure unique est une representation 
schematique d'une mise en oeuvre preferee de la methode selon 
r invention. 

Le forage d'un puits 1 permet d'accumuler un grand nombre 
d' informations sur la zone du sous-sol dans lequel est for6 le puits. Des 
donnees sismiques sont 6galement enregistrees independamment et dans la 
meme zone, et sont utilis6es notamment pour la constmction d'un bloc 
somme 2 appel6 dans ce qui suit "bloc somme enregistre". 

Le bloc somme enregistre 2 est utilise pour Teiaboration d'un bloc 
de donnees sismiques migrees-temps' 3, generalement k partir d'un 
modele de vitesses concemant la zone consideree, appeie modele ou 
champ de vitesses de migration-temps 13, represente en fonction d'une 
echelle verticale appcjee 'temps migre'. Ce moddle est une representation 
trfes grossidre du sous-sol, qui peut etre represente par son champ de 
vitesses de propagation des ondes sismiques V(x,y,z) incluant les 
heterogeneites ainsi que I'anisotropie symbolisee par le champ E(x,y,z), 
mal connus du geophysicien. Toute amelioration de la connaissance de 
cette information sur la realite geplogique du sous-sol est un gain 
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important pour T exploration. Ces techniques sont bien connues des 
specialistes et ne seront pas d6crites dans le detail. 

Les informations recueillies dans Ic forage 1 doivent 6tre 
confrontees, dans la pratique courante d' exploration, aux donn£es 
sismiques enregistrdes. 

Cette confrontation se fait generalement dans le domaine des 
donnees sismiques migr6es-temps. 

Le plus souvent, on construit en partant des donnees sismiques 
migr6es-temps ct de toutes les informations relatives de pres ou de loin 
aux vitesses de propagation (ten^s de trajet, vitesse de sommation, etc.) 
un module vitesse-profondeur 4, geologiquement representatif de la zone 
consid6r6e, de nombreuses techniques plus ou moins performantes mais 
toutes approximatives ctant disponibles pour r6aliser ce qui est 
communement appele la conversion profondeur ou la migration 
profondeur selon la technique de mise en oeuvre. 

En eflfet, les informations de puits (position mais surtout pendage) 
sont connues avec une precision limit6e, les vitesses sont biaisees par le 
phenomene d'anisotropie et les proced6s d'daboration des modeles de 
vitesses sont insuffisamment precises. 

La presente invention a pour objet principal une methode pour caler 
la trajectoire du puits, repr6sent6e dans le moddle de vitesse-profondeur 4, 
dans le bloc sismique migr6-temps, ce qui n'a pu 6tre r6alis6 jusqu'il 
pr6sent avec pr^ision, en raison justement de toutes les hypothdses ou 
approximations qui sont faites dans les techniques disponibles. 

Quelle que soit la technique utilis6c pour T elaboration du modde 
de vitesse-profondeur 4, comme par exemple un 6tirement vertical d'un 
bloc migre-temps, il est necessaire de caler ledit module de vitesse- 
profondeur 4 sur le puits 1 et de lui associer la trajectoire 5 du puits dans 
le modele vitesse-profondeur 4. Ainsi chaque point A de la trajectoire 5 
est d6fini au moins par des coordonnees spatiales X, Y, Z et un vecteur 
pendage (pendage et azimut) 

Le vecteur pendage pent 6tre obtenu directement a partir du modele 
vitesse-profondeur 4 lorsque le point A est a T intersection d'un r6flecteur 
et de la trajectoire. Mais, ce vecteur pendage pent 6tre obtenu par diverses 
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voies telles que la pendagemetric dans le puits, T interpolation ou 
extrapolation entre deux rdflecteurs profondeurs connus. 

Le modele vitesse-profondeur 4 est repr6sent6 sur la figure sous la 
forme d'un modele k couches definies par leur champ de vitesse Vf, V2, 
5 V3 et qui sont d61imit6es entre elles par des rdflecteurs dont trois 6^8 
sont representds. Mais il pourrait etre repr^sente sous la forme maill6e, 
representation qui est bien connue encore des sp^cialistes, la forme 
maillee etant obtenue a partir d*une grille dans le sous-sol et dont les 
noeuds affect6s d'une valeur de vitesse presentent un int6ret pour 
10 I'elaboration du reseau maille dans le modele vitesse-profondeur 4. 

Dans le modele vitesse-profondeur 4, la U-ajectoire 5 intersecte les 
reflecteurs 6 a 8 appel6s reflecteurs profondeur. Chaque r£fi[ecteur 
profondeur 6 a 8 est d6fini enure autre par son pendage et son azimut. 

Dans une mise en oeuvre pr6feree de Tinvention, on s'int6resse 
15 plus particulierement aux points d' intersection entre la trajectoire 5 et les 
reflecteurs profondeurs 6 a 8 et la methode va etre decrite pour Tun de ces 
points d' intersection, par exemple ie point A, 6tant precise que la ra6thode 
peut etre appliqu6e a tous les auu-es points de la trajectoire et aux points 
situes autour de chacun de ces points d* intersection reels ou presumes, 
20 comme on le verra ultdrieurement. 

Pour chaque point A du modele vitesse-profondeur 4, on procdde k 
une mod^lisation par lancer du rayon k offset ou deport nul, k travers ledit 
moddle de vitesse-profondeur 4 ; le rfisultat de cettc operation est un point 
image a, qui est dafini par ses coordonnSes spatio-ten^)orelles et par un 
25 couple pendage-temps dans un domaine synthetique 9 equivalent au bloc 
somme enregistre 2. 

Pour diverses valeurs de pendage et d'azimuts du point A, le lancer 
de rayons a offset nul permet de construire, dans le domaine synthetique 
9, une surface de diffraction 10 sur laquelle sont situ6es toutes les images 
30 a du point A, pour les valeurs de pendage et d'azimut considerees. 

Dans une autre 6tape, on superpose le domaine synthetique 9 au 
bloc somme cnregistr6 2, cette operation de superposition etant 
materialisee par une fldche 1 1 sur la figure. 

Le resultat de la superposition est une possible coincidence ou 
35 tangence de la surface de diffraction 10, representee en pointilie dans le 
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bloc somme enregistre 2, avec Ic r6flecteur 12. Le point de tangence ou 
de coincidence entre la surface de diffraction 10 et le r6flecteiir 12 est 
appeia point spaculaire a. Comme il y a unc relation biunivoque entre le 
point speculaire a et le pendage correspondant au point A, ce pendage est 
une information nouvelle sur le pendage r6el du r6flecteur au point A. 

II est egalement possible de completer la surface de diffraction 10 
pour tenir compte par exemple de T incertitude de mesure de la trajectoire 
spatiale du puits, en prenant plusieurs points ou un nuage de points situes 
autour et proches du point d' intersection seIectionn6 A et de proceder a 
des lancers de rayons a offset nul pour chacun de ces points afm de mieux 
pr^ciser la position du point sp6culaire a. 

La meme operation est repet^e pour chacun des points 
d' intersection ou proches de T intersection de la trajectoire 5 avec les 
reffecteurs profondeurs 6^8, ...etc. du modele 4, chacun des points 
d' intersection A 6tant bien evidemment rep6r6 par ses cooidonn6es X^, 
Ya, Za- Lorsque le point A n'appartient pas a un r6flecteur profondeur, 
on prend des points situes de part et d'autre de la trajectoire 5 et autour du 
point A et par interpolation on d6termine les parametres affectes au point 
A, puis on lui applique T operation pour determiner son image a dans le 
domaine synth6tique 9. 

A cette etape, il est possible de proceder a une verification visuelle 
et k une correction 6ventuelle du modele vitesse-profondeur 4 et ce, de 
manidre interactive ou en temps rfiel pendant la mise en oeuvre de la 
methode. 

En effet, il devient possible de procfider sur la station de travail k 
une comparaison entre le domaine somme synth6tique 9 et le bloc somme 
enregistr6 2 de maniere a determiner si un ecart a la tangence ou i la 
coincidence existe entre chaque surface de diffraction 10 et le r6flecteur 
correspondant 12 du bloc somme enregistre. Lx)rsque cet ecart existe, on 
procede a des iterations successives sur le modele vitesse-profondeur 4, 
par changement des valeurs du vecteur pendage-azimut et/ou des valeurs 
de vitesses du module vitesse-profondeur et/ou des paramdtres 
caract6risant Tanisotropie de la zone du puits, jusqu'i ce que ledit dcart 
soit acceptable pour considerer que ledit modele vitesse-profondeur 4 est 
correctement cal6. Ainsi, il est possible d'afFmer le moddle vitesse- 
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profondeur 4 au cours de la mise en oeuvre de la in6thode selon 
r invention et de caler avec plus de precision la trajectoire 5 du puits dans 
ledit modele. 

Cela pennet egalement de verifier si rintervalie de pendage qui a 
ete utilise est correct pour le point A selectionn6. 

Dans une 6tape suivante, on procdde a une migration temps pour 
chacun des points de chaque surface de diffraction 10 du domaine somme 
synth6tique 9, par exemple par simulation, en tenant compte du couple 
pendage- temps correspondant a chacun de ces dits points, laquelle 
simulation est' effectuee k travers le champ de vitesse ten^)s-migre 13. 
Pour realiser cette simulation au moins deux proc6d6s peuvent dtre 
appliques. 

Un des proc6des est celui mettant en oeuvre un lancer de rayons 
anisotropes dans le modele de vitesses de migration temps (B. 
RAYNAUD 1993 : Real time migration operators simulated by 
anisotropic ray tracing. 55th EAEG-C045), pas a pas, pour simuler 
r operation de migration temps par differences finies ou tout autre procdde 
prenant en compte de maniere 6quivalente le modele de vitesses de 
migration temps. 

Un autre proc6d6 consiste a simuler la migration temps de 
KIRCHHOFF par une approche analytique. Ce proc6d6 consiste k migrer 
un point de la surface de diffraction suppos6e Stre un hyperboloide pour 
I'amener a I'apex dudit hyperboloide. 

Dans ces procdd6s de simulation de la cinematique de la migration- 
temps, on peut utiliser avantageusement une loi de vitesse qui tient compte 
des variations de vitesse incluses dans le modele de vitesse de migration- 
temps de la zone consideree. 

Le resultat de la simulation dans un domaine migration-teii^}S 14 de 
chacune des surfaces de diffraction 10 n'est pas un point unique mais 
encore une surface migree-temps 15, d6nomm6e ici "plumei" (plume en 
Anglais). 

Comme on Ta ddcrit pour le domaine sonmiation, on peut 
Egalement dans le domaine temps-migre superposer le domaine 14 au bloc 
temps-migr6 3, cette operation de superposition 6tant mat6rialis6e par une 
fidche 16 sur la figure. 
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Le r6sultat de la superposition est une possible coincidence ou 
tangence du plumet 15 repr6sent6 en pointili6 dans le bloc migr^-temps 3 
qui conqirend des horizons sismiques niigr6s-temps 17, 18. L* image du 
point a est en a' dans le domaine 14 et celle du point sp6culaire est en a' 
dans le bloc 3. 

A cette etape, il est encore possible de proceder a une verification 
visuelle et a une eventuelle correction du modele vitesse-profondeur 4 et 
ce, de manidre encore interactive ou en temps r6el pendant la mise en 
oeuvre de la m6thode. En effet, on peut comparer chaque plumet 15 
migr6-temps et un 6ventuel horizon sismique 18, de manidre k d6terminer 
s'il y a un 6cart a la tangence ou coincidmce entre lesdits plumet et 
horizon sismique. Lorsque cet ccart existe, on procede k des iterations 
successives sur le modele vitesse-profondeur 4, par changement des 
valeurs du vecteur pendage-azimut et/ou des valeurs de vitesse du modele 
vitesse-profondeur et/ou des parametres caract6risant T anisotropic de la 
zone du puits, jusqu'a ce que ledit ecart soit acceptable pour considerer 
que ledit modele vitesse-profondeur 4 est correctement cal6. 

Ainsi, la presente invention penriet sensiblement k chaque 6tape de 
la m6thode de procMer a des corrections des d^fauts 6ventuels, de revenir 
aux causes desdits defauts et d'y rem6dier par des actions ponctuelles et 
ce, tras rapidement, le rasultat final 6tant un excellent et pracis calage de 
la trajectoire 5 dans le bloc migra-temps et une amaiioration du moddle 
migration apres et avant sommation des donnaes sismiques, imageries, 
ten^)s-profondeur, etc. et du modele vitesse-profondeur 4, ce modele 
devenant alors apte a etre utilise a d'autres fins. 

Dans rhypothese ou la connaissance du champ de vitesse 13 serait 
imparfaite, il est possible selon T invention d'ameliorer ladite 
connaissance. En effet, k partir d*un pointa dans le bloc migra-temps 3, 
par exemple des horizons 17, 18, on crae des raflecteurs correspondants 
qui peuvent etre transformas en des raflecteurs somme par des processus 
de damigration strictement inverses des simulations de migration-temps. 
Ces raflecteurs damigras sont alors superposas aux raflecteurs ou donnaes 
sismiques du bloc 2. Tout acart est nota et peut etre utilisa pour amaiiorer 
le champ de vitesse 13 . 
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Une autre caract6ristique de la pr6sente invention est de pennettre 
une mise en oeuvre inverse de la m6thode. En effet, pour un point de 
tangence entre un plumet 15 et Thorizon correspondant dans le bloc 
migre-temps, on peut d6tenniner le couple pendage-temps qui lui est 
associ6. Puis, on repere dans le domaine somme synth6tique 9 la surface 
de diffraction 10 correspondant audit plumet 15 et son point sp^culaire a 
dans le bloc sommation 2. Comme le point speculaire a correspond k un 
seul point du modele vitesse-profondeur 4, il devient possible de 
d6ternuner la position de ce point ainsi que le vecteur pendage-azimut 
associ^. 

Grice S la pr&ente invention, on obtient, pour un module 
profondeur 4 choisi et pour un point particulier A de la trajectoire 5 et 
pour un vecteur-pendage s61ectionne, les positions correspondantes a du 
bloc sommation 2 et a' dans le bloc migre-temps 3. 

S'il y a ambiguite ou incertitude sur le pendage, le point A se 
projette en une autre surface a" (non representee) caract^ristique de cette 
ambiguite. 

La superposition dans Tun ou I* autre domaine sonunation ou 
migration-temps perment 6venttiellement de lever cette ambiguity, mais 
pas obligatoirement. Dans ie dernier cas le resultat du proc6de n'est pas 
une trajectoire unique, mais une "zone-trajectoire" possible. 

Un avantage de la pr6sente invention qui decoule de ce qui pr6c&le 
est que soit a partir d*un point ou d*une zone de la surface de dif&action 
10, soit k partir d'un point ou d'une partie du plumet 15, on peut revenir 
de fa(on constante au moddle de vitesse-profondeur 4. 

II est ais6 de comprendre egalement qu'im point de tangence ou de 
coincidence sur les blocs 2 et/ou 3, permet de revenir directement au 
moddle vitesse-profondeur 4 sans etape intermediaire et d'en ddduire des 
informations trds pr6cises et tres utiles aux sp6cialistes. 
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REVENDICATIQNS 

1. Methode pour representer dans un domaine spatio-temporel la 
trajectoire (5) d'au moins un puits (1) fore dans une zone du sous-sol, 
caracteris6e en ce qu*elle consiste a 

a) 61aborer un modele vitesse-profondeur (4) cale au puits (1), 
geologiquement reprdsentatif de ladite zone et comportant une 
representation de iadite trajectoire (5) dudit puits ; 

b) seiectionner au moins un point (A) sur ou pr&s de ladite 
trajectoire, ledit point (A) 6tant d6fini par ses coordonnees spatiales et par 
au moins un couple de donn^es pendage-azimut ; et 

c) au moins modeliser dans un domaine sonune synthetique (9) 
rimage (a) du point selectionne (A), ladite image (a) 6tant d'une part 
definie par des coordonnees spatio-temporeiles et par un vecteur pendage- 
temps et, d' autre part, obtenue par un lancer de rayon a deport nul dans le 
moddle vitesse-profondeur (4). 

2. Methode selon la revendication 1, caracterisee en ce qu*en outre 
elle comprend Tetape d'effectuer une migration en temps de ladite image 
(a) du domaine synthetique (9), de maniere k obtenir une representation 
(a') de ladite image (a) dans un domaine migre-tenq)s (14) 

3. Methode selon la revendication 1, caracterisee en ce qu'en outre le 
domaine sonmie synthetique (9) contenant Timage (a) est superpose k un 
bloc somme (2) de donnees sismiques. 

4. Methode selon la revendication 2, caracterisee en ce qu'en outre le 
domaine migre-temps (14) contenant la representation (a') est superpose i 
un bloc de donnees sismiques migrees-temps (3). 

5. Methode selon la revendication 4, caracterisee en ce que le bloc de 
donnees sismiques migrees-temps est obtenu a partir du bloc somme 
enregistre (2) en utilisant un modele de vitesse migration donne (13) et en 
ce que la migration de ladite image (a) est effectuee au moyen dudit 
moddle de vitesse (13). 

6. Methode selon la revendication 1, caracterisee en ce que module de 
vitesse-profondeur (4) est a couches deiimitees par des r6flecteurs (6 k 8). 

7. Methode selon la revendication 1, caracterisee en ce que le modele 
de vitesse-profondeur (4) est un modele mailie. 
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8. M6thode selon les revendications 1 et 6, caract6ris6e en ce qu'on 
seiectionne plusieurs points (A) dans le module vitesse-profondeur (4) et 
plus particulidrement les points situes i ou proches de T intersection de la 
trajectoire (5) avec les r6flecteurs (6 k 8). 

9. M6thode selon Tune des revendications 1 k 8, caract6ris6e en ce 
que le bloc de donnees sismiques migrees-temps (3) est obtain par 
migration de KIRCHHOFF du bloc de donndes sismiques somme (2) et en 
ce que le domaine migra-temps (14) resulte d'une sommation analytique 
de ladite migration de KIRCHHOFF. 

10. Mediode selon Tune des revendications 1 a 8, caracterisee en ce 
que bloc de donnees sismiques migr£es-temps (3) est obtenu au moyen 
d'une migration par differences finies et en ce que le domaine migre- 
temps rasulte d'une operation de lancer de rayons anisotropes pas pas 
dans le domaine de vitesse (13). 

1 1 . Methode selon Tune des revendications precddentes, caract6ris6e en 
ce que chaque image (a) dans le domaine somme synthetique (9) est 
distribute sur une surface dite surface de diffraction (10). 

12. Mtthode selon la revendication 11, caracterisee en ce qu'on 
s61ectionne une partie de la surface de diffraction (10) comme 6tant 
rimage dans le domaine somme synthdUque (9) du point (A) selectionn6 
dans le module de vitesse-profondeur (4). 

13. M6thode selon Tune des revendications pr6c&lentes, caract6ris6e en 
ce qu'on rtalise, au cours de la mise en oeuvre de ladite mfithode, une 
comparaison entre le domaine somme synthetique (9) et le bloc somme 
enregistre (2) et on determine Tteart a la tangence ou k la coincidence 
susceptible d'apparaitre entre chaque surface de diffraction (10) et le 
reflecteur (12) correspondant et, lorsqu'il y a un ecart, on procede k des 
iterations successives sur le modele vitesse-profondeur (4) jusqu'a ce que 
ledit 6cart soit acceptable pour considerer que ledit modele vitesse- 
profondeur (4) est correctement cal6. 

14. Methode selon la revendication 11, caracttris^e en ce que la 
representation en temps migrfi de la surface de diffraction (10) est une 
surface migr6e-temps (15). 

15. Methode selon la revendication 14, caracterisee en ce qu'on 
selectionne une partie de la surface migree-tcn^>s (15) comme 6tant 
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r image dans le domaine migre-temps (14) du point (A) s6Iectionne dans ie 
moddle vitesse-profondeur (4). 

16. M6thode selon la revendication 14, caract6ris6e en ce qu'au cours 
de la mise en oeuvre de ladite methode, on d6termine s*il y a un 6cart k la 

5 tangence ou k la coincidence entre ladite surface migree-temps (15) et 
r horizon migre-temps correspondant (18) et, lorsqu'il y a un 6cart, on 
procdde a des iterations successives sur le modfele vitesse-profondeur (4) 
jusqu'a ce que ledit 6cart soit acceptable pour considerer que ledit modele 
vitesse-profondeur (4) est correcteraent cal6. 

10 17. Methode scion Tunc des revendications 13 ou 16, caract6risee en ce 
que le point de tangence (a ou a'), lorsqu'il existe, permet de d6tenniner 
le pendage rtel du r6flecteur dans le module vitesse-profondeur (4). 
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